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Abstract 



We study the decay of scalar fields, in particular the inflaton, into lighter 
scalars in a De Sitter spacetime background. After providing a practical def- 
inition of the rate, we focus on the case of an inflaton interacting with a 
massless scalar field either minimally or conformally coupled to the curvature. 
The evolution equation for the expectation value of the inflaton is obtained to 
one loop order in perturbation theory and the decay rate is recognized from 
the solution. We find the remarkable result that this decay rate displays an 
equilibrium Bose-enhancement factor with an effective temperature given by 
the Hawking temperature H/2w, where H is the Hubble constant. This con- 
tribution is interpreted as the "stimulated emission" of bosons in a thermal 
bath at the Hawking temperature. In the context of new inflation scenarios, 
we show that inflaton decay into conformally coupled massless fields slows 
down the rolling of the expectation value. Decay into Goldstone bosons is 
also studied. Contact with stochastic inflation is established by deriving the 
Langevin equation for the coarse-grained expectation value of the inflaton field 
to one-loop order in this model. We find that the noise is gaussian and cor- 
related (colored) and its correlations are related to the dissipative ("decay") 
kernel via a generalized fluctuation-dissipation relation. 
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I. INTRODUCTION 



Non-equilibrium processes play a fundamental role in cosmological scenarios motivated 
by particle physics models. Processes such as thermalization, reheating and particle decay 
are important ingredients in most theoretical scenarios that attempt a description of early 
universe cosmology An important ingredient underlying most attempts at a descrip- 

tion of inflationary cosmologies is that of equilibration, which is tacitly assumed whenever 
quasi-equilibrium methods are invoked in their study, such as effective potentials and finite 
temperature field theory. 

The concept of quasi-equilibrium during an inflationary stage requires a careful under- 
standing of two different time scales. One is the expansion time scale tjg — H^ 1 , where 
H is the Hubble parameter, while the other is the field relaxation time scale, t r . A quasi- 
equilibrium situation obtains whenever t r «tu since in this case particle physics processes 
occur much faster than the Universe expands allowing fields to respond quickly to any changes 
in the thermodynamic variables. 

Can the inflaton of inflationary cosmologies be treated as being in quasi-equilibrium? 
The answer to this question depends crucially on the various processes that contribute to 
the dynamical relaxation of this field. Some of these processes are the decay of the inflaton 
into lighter scalars or fermions, as well as collisions that, while they leave the total number 
of particles fixed, will change the phase space distribution of these particles. 

In this paper we wish to understand the effects of coupling the inflaton to lighter 
scalars during the inflationary era. Single particle decay of the inflaton into lighter scalars 
or fermions in the post-inflationary era, when space-time can be well approximated by a 
Minkowski metric, is by now fairly well understood within the context of the "old-reheating" 
scenario |6|-|9|], but it is much less understood during the early stages of inflation. 

In saying this, we have to distinguish between the various models of inflation. In "old" or 
"extended" inflation |H| , the inflaton is trapped in a metastable minimimum, from which it 
exits eventually via the nucleation of bubbles of the true vacuum phase. A critical parameter 
in these models is the nucleation rate, which is typically calculated assuming that the field 



is equilibrated p|,^|Jll|. Such an assumption will only be justified if the relaxation time of 
the inflaton in the metastable phase, i.e. in De Sitter space, is much smaller than the time 
scale for nucleation of a critical bubble. 

In "new" or "chaotic" inflationary scenarios |7|], the supraluminal expansion of the uni- 
verse is driven by the so-called "slow-roll" dynamics of the inflaton. Typically, this is studied 
via the zero or finite temperature effective potential. In doing this, the assumption of quasi- 
equilibrium is made, in that one can replace the actual field dynamics, by that given by 
a static quantity, the effective potential. This assumption must be checked, and one way 
of doing that is by considering the evolution of the inflaton in the presence of couplings to 
other fields and computing the field relaxation rate. In particular, if dissipative processes 
associated with the "decay" of the inflaton field into lighter scalars are efficient during the 
slow-roll stage in new or chaotic inflationary scenarios, then these may result in a further 
slowing down of the inflaton field as it rolls down the potential hill. 
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Finally, there is the stochastic inflation approach advocated by Starobinskii ||12||. In this 



scenario, inflationary dynamics is studied with a stochastic evolution equation for the inflaton 
field with a gaussian white noise that describes the microscopic fluctuations. As pointed out 
by Hu and collaborators [13| and Habib and Kandrup [14|, IS], noise and dissipation are 



related via a fundamental generalized fluctuation-dissipation theorem and must be treated 
on the same footing. 

In this article we study the process of relaxation of the inflaton field via the decay into 
massless scalar particles in De Sitter space-time with the goal of understanding the modifi- 
cations in the decay rates and the relevant time scales as compared to those in Minkowski 
space-time. Our study focuses on the description of these processes within "old/extended 
" and "new" inflationary scenarios. Within "old" inflation we study the situation in which 
the inflaton field oscillates around the metastable minimum (before it eventually tunnels) 
and decays into massless scalars. In the case of "new" inflation, we address the situation of 
the inflaton rolling down a potential hill, near the maximum of the potential, to understand 
how the process of "decay" introduces dissipative contributions to the inflaton evolution. 

We also establish contact with stochastic inflation by deriving the Langevin (stochastic) 
equation for the inflaton field, and analyzing the generalized fluctuation-dissipation relation 
between the dissipation kernel resulting from the "decay" into massless scalars and the 
correlations of the stochastic noise. 

At this point we would like to make precise the meaning of "decay" in a time dependent 
background. The usual concept of the transition rate in Minkowski space, i.e. the transition 
probability per unit volume per unit time for a particular (in) state far in the past to 
an (out) state far into the future is no longer applicable in a time dependent background. 
Relaxing this definition to finite time intervals ( for example by obtaining the transition 
probability from some t« to some t/) would yield a time dependent quantity affected by the 
ambiguities of the definition of particle states. In Minkowski space-time the existence of 
time translational invariance allows a spectral representation of two-point functions and the 
transition rate (or decay rate) is related to the imaginary part of the self-energy on shell. 
Such a correspondence is not available in a time dependent background. Furthermore the 
time evolution of a scalar field is damped because of the red-shifting of the energy due to the 
expanding background even in the absence of interactions. In particular in De Sitter space 
this damping is exponential and any other exponential damping arising from interactions 
will be mixed with the gravitational red-shift. 

We propose to define the "decay rate" as the contribution to the exponential damping of 
the inflaton evolution due to the interactions with other fields, which is clearly recognized 
through the dependence on the coupling to these fields. This practical definition then argues 
that in order to recognize such a "decay rate" we must study the real-time dynamics of the 
non-equilibrium expectation value of the inflaton field in interaction with other fields, and 
the identification of this rate should transpire from the actual evolution of this expectation 
value. 

Although the methods presented in this article can be generalized to arbitrary spatially 
flat FRW cosmologies and include decay of a heavy inflaton into lighter scalars, we study 
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the simpler but relevant case of inflaton decay into massless scalars in De Sitter space-time. 
This case affords an analytic perturbative treatment that illuminates the relevant physical 
features. 

The article is organized as follows: in section II we present a brief pedagogical review of 
the closed time path formalism and the derivation of the non-equilibrium Green's functions 
for free fields in De Sitter spacetime. In section III the model is defined and the methods 
to obtain the evolution equation are summarized. Here we study several cases, including 
minimal and conformal coupling and the slow-roll regime. We obtain the remarkable result 
that the decay rate can be understood from stimulated emission in a bath at the Hawking 
temperature. In section IV we study the effect of Goldstone bosons on the inflaton evolution, 
pointing out that Goldstone bosons are necessarily always minimally coupled. In section V 
we make contact with stochastic inflation, derive the stochastic Langevin equation for the 
inflaton zero mode and obtain the dissipative kernel and the noise correlation function, 
which are related by a generalized fluctuation-dissipation theorem. The resulting noise is 
gaussian but colored. In section VI we summarize our results and propose new lines of 
study. An appendix is devoted to a detailed derivation, establishing the consistency of the 
flat-space limits of the analytic structure of the dissipative kernel in De Sitter space-time, 
with Minkowskian calculations. 



II. NON-EQUILIBRIUM FIELD THEORY 

A. The Closed Time Path Formalism 

As mentioned in the introduction we will be studying the real-time evolution of the zero 
mode of the inflaton in a universe with a time-dependent background metric. Furthermore, 
we will restrict our investigations to perturbative phenomena in theories in which the inflaton 
interacts with other matter fields in a De Sitter universe. We will be primarily concerned with 
defining and calculating quantities such as the perturbative "decay rate" of the inflaton. For 
studies of non-perturbative particle production in non-equilibrium situations see |19| (and 
references therein). Calculation of particle production in curved spacetime has also been 
done using the Boguliubov transformations; see for e.g | p3 ] 



The tools required for studying time-dependent phenomena in quantum field theories have 
been available for quite some time, and were first used by Schwinger Jl6| and Keldysh ||17|| . 



There are many articles in the literature using these techniques to study time-dependent 
problems pT8|JT9|] . However, these techniques have not yet become an integral part of the 



available methods for studying field theory in extreme environments. We thus present a 
concise pedagogical introduction to the subject for the non-practitioner. 

The non-equilibrium or time dependent description of a system is determined by the time 
evolution of the density matrix that describes it. This is in turn described by the Liouville 
equation (in the Schrodinger picture): 

ih^ = [H(t),p(t)\. (2.1) 
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Here we have allowed for an explicitly time dependent Hamiltonian which is in fact the case 
in an expanding universe. The formal solution of the Liouville equation is, 

p{t) = U{t) Pl U-\t) (2.2) 

where pi is the initial density matrix specified at some initial time t = 0. The time dependent 
expectation value of an operator is then easily seen to be, 

(0)(t) = TrfaU-WOUWyTrlpi] (2.3) 

In most cases of interest the initial density matrix is either a thermal one or a pure state 
corresponding to the ground state of some initial Hamiltonian. In either case, 

Pl = e-P Hi (2.4) 
Hi = H(t < 0). (2.5) 

The ground state of Hi can be projected out by taking the j3 — ► oo limit. To study strongly 
out of equilibrium situations it is usually convenient to introduce a time-dependent Hamilto- 
nian H(t) such that H(t) = Hi for — oo < t < and H(t) = H evo i(t) for t > where H evo i(t) 
is the evolution hamiltonian that determines the dynamics of the theory. 

For thermal initial conditions, eq ( |2.3| ) can be written in a more illuminating form by 
choosing an arbitrary time T < for which U(T) = exp[— iTHj\ so that we may write 
exp[-(3Hi] = exp[-iHi(T-ip-T)] = U(T-i/3,T). Inserting the identity U-\T)U(T) = 1, 
commuting U~ l {T) with pi and using the composition property of the evolution operator, 
we get 

(0)(t)= Tr[U(T - 0, t)OU(t, T)}/Tr[U(T - i/3, T)] (2.6) 



= Tr[U(T - i/3, T)U(T, T')U{T', t)OU(t, T)}/Tr[U(T - i/3, T)] (2.7) 

where we have introduced an arbitrary large positive time T' . The numerator now represents 
the process of evolving from T < to t, inserting the operator O, evolving to a large positive 
time T', and backwards from T' to T and down the imaginary axis to T — i(3i. Eventually 
one takes T — ► — oo; T' — > oo. This method can be easily generalised to obtain real time 
correlation functions of a string of operators. 

As usual, an operator insertion may be achieved by introducing sources into the time 
evolution operators, and taking a functional derivative with respect to the source. Thus we 
are led to the following generating functional, in terms of time evolution operators in the 
presence of sources 

Z[J + , J~, J 13 } = Tr[U(T - i(3, T; J P )U(T, T'; J~)U(T', T; J + )] (2.8) 

with T — ► — oo;T' — > oo. One can obtain a path integral representation for the generating 
functional [|I6H 21~H by inserting a complete set of field eigenstates between the time evolution 
operators, 
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Z[J+, J~, .f] = J DED^D^ J VZ+VE-VEPjfi {m + ,J + ]~m ,j ]} x ^ ^ 

The above expression for the generating functional is a path integral defined on a complex 
time contour. Since this path integral actually represents a trace, the field operators at the 
boundary points have to satisfy the following conditions: 

E+(T) = Tf(T-ip) = £, (2.10) 
E+(T') =£-(T') = £ 2 , 
£~(T) = £^(T) = Ei. 

The superscripts (+) and (— ) refer to the forward and backward parts of the time contour, 
while the superscript {(3) refers to the imaginary time part of the contour. In the limit 
T — > — 00 it can be shown that cross-correlations between fields defined at real times and 
those defined on the imaginary time contour, vanish due to the Riemann — Lebesgue lemma. 
The fact that the density matrix is thermal for t < shows up in the boundary condition 
on the fields and the Green's functions. Thus the generating functional for calculating real 
time correlation functions simplifies to: 

Z[J + , J~\ = exp \i£ dt[C int {-i5/5J + ) - C int (i5/5J-)]\j x 

dhJ T dt 2 J a {t x )J b {t 2 )G ab {t x ,t 2 )\ (2-11) 

with a, b = +, — . We will now proceed to a calculation of the non-equilibrium Green's 
functions in De Sitter spacetime 




B. Non-equilibrium Green's Functions In De Sitter 

To obtain the non-equilibrium Green's functions we focus on a free scalar field £ in 



an FRW background metric [20]. The generating functional of non-equilibrium Green's 
functions is written in terms of a path integral along the closed time path (CTP) ( j2.11j) with 
a free Lagrangian density, £o(S ± ) given by 

£ (E ± ) = i[a 3 (t)E+ 2 - a(t)(VE+) 2 - a 3 (t)(M 2 + ^)S +2 ] (2.12) 



where 1Z = 6(a/a + a 2 /a 2 ). 

We prepare the system such that it is thermal for t < with temperature 1/(3. At t — 0, 
we switch the interactions on and follow the resulting time evolution of the coupled fields. 
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As seen in the previous section, the CTP formulation of non-equilibrium field theory imposes 
certain boundary conditions on the fields ( |2.10|) . 

Rewriting the Lagrangian as a quadratic form in E we obtain the Green's function 
equation, 

§ + 41 " ? + < M2 + ^ G <*> i; W = M0¥Y (2 - 13) 

We use spatial translational invariance to define, 

G k {t,t') = Jd 3 xe ils G{x,t-0,t'), (2.14) 
so that the spatial Fourier transform of the Green's function obeys, 

£ + 3 - a j t + 1 + (M2 + * - sfefe- (2 - 15) 

The solution can be cast in a more familiar form by writing 

GkM = ^$^kt r y (2,16) 

so that the function f k (t,t') obeys the second order differential equation, 
The general solution of this equation is of the form 

f k (t,t') = f k >(t,t')e(t - f) + /*<(*, oe^ - 1), (2.18) 

where and f k K are solutions to the homogeneous equation obeying the appropriate 
boundary conditions. They can be expanded in terms of normal mode solutions to the 
homogeneous equation: 

/*>(*, = A>U h (t)U*(t') + B>U k (t')U*(t), (2.19) 
/*<(*, = A<U k (t)U* k (t') + B<U k (t')U* k (t). 

The choice of mode functions can be made unique for a given £ and M by requiring that, 

-ikt 

-ii u k yL) ~ e 

»oo 

or 



lim U k (t) ~ e~ iH , (2.20) 



hm £7 fc (t) ~ e 



u k = Vk 2 + M 2 . 
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Here H = a(t)/a(t) is the Hubble parameter. Stated simply, these conditions require that we 
must necessarily recover the rules of field theory in Minkowski space when we take the flat- 
space or the short distance (k — > oo) limits. These boundary conditions are similar to those 
invoked by Bunch and Davies |22|.|2l|]. For the De Sitter space case, with a(t) = exp(Ht), 



the mode functions are the Bessel functions [p3| , J± u (ke / H) provided v is not an integer, 
where 



M 2 9 

- = V^ + 12 ^4- (2 ' 21) 

For the values of £ and M for which v becomes real these Bessel functions become manifestly 
real and the solutions that fulfill the boundary conditions (|2.20|) are obtained as linear 
combinations of these functions. In order to obtain the full Green's functions the constants 
A > ,B > ,A < ,B < ( |2.19 ) are determined by implementing the following boundary conditions 



f>(t,t) = f<(t,t) (2.22) 

f>(t,t)-f<(t,t) = l 

f>(T-i(3,t)=f<(T,t) 

The first two conditions are obtained from the continuity of the Green's function and the 
jump discontinuity in its first derivative. The third boundary condition is just a consequence 
of the periodicity of the fields in imaginary time, which followed from the assumption that 
the density matrix is that of a system in thermal equilibrium at time T < 0. Although 
we allowed for a thermal initial density matrix, we will focus only on the zero temperature 
case, since during the De Sitter era the temperature is rapidly red-shifted to zero. The most 
general case for the Green's function can be found in ||20|| . Implementing these conditions 
for the case of imaginary v and zero temperature (j3 — > oo) we find : 

r>(+ +'\ — -^Ju{z)J-y{z') . . 

k 1 ' } ~ 2tfsinMe 3 ^/2 e 3ffi'/2 [Z - Z6) 



2#sin(7rz/)e 3ffi /2 e 3Hi'/2' 



z = fc— , z' = k— (2.24) 

and v is given by ( |2.21| ). 

In the case of interest for this article, that of a massless scalar field either conformally or 
minimally coupled to the curvature, the Bessel functions in the above expressions must be 
replaced by the corresponding Hankel functions. These two cases correspond to: 

[i)M = 0, £ = 0, U k (t) = H^ 3/2 (ke- Ht /H), (2.25) 

{ii)M = 0, £ = ~ U k (t) = H^ l/2 (ke- m /H), (2.26) 
o 
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Finally, the Green's functions for these two cases of a massless scalar field, in terms of z and 
z', are: 

-inHU(z)H<?)(2>) 

C<(t t') - -iKHPWWW 

with v — 3/2 for minimally coupled or v — 1/2 for conformally coupled. From the expressions 
for G^T and we can compute the required 2-point functions: 

i(TS+(f, i)£+(r , t')) = G > {f, t; 7, t')G(t - t') + G < (f, t; 7, t')Q(t' - t), (2.29) 
i(TE-(r, t)S~(r , t')) = G > (f, t; f 1 , t')Q(t' - t) + G < (f, t; 7, t')G(t - if), (2.30) 
-i(TE + ^(f,0S- (+) (^,O) = G <{>) (r,t;7,t'). (2.31) 

In our subsequent analysis we will always work in the f3 — > oo limit (zero initial temper- 
ature) unless stated otherwise. 



III. THE MODEL AND THE METHOD 

We now turn to the main topic of this work, the study of an inflaton, described by a field 
$ coupled to a massless scalar field a with the non-equilibrium Lagrangian density 

C = I[e 3 ^$+ 2 - e ^(V<l» + ) 2 - e mt (ml + 12£,# 2 )$ + 2 (3.1) 
+e ^t &+ 2 _ e m^ a +f _ i2e™ t i a H 2 a +2 ] + ^e 3Ht [g^ + a +2 + h(t)$ + ] 

-[+ — -]■ (3-2) 
We have introduced a "magnetic" (source) term h(t) for two reasons: 

1. such a term will be generated by renormalization from tadpole graphs, therefore we 
introduce h(t) as a counterterm to cancel these contributions, 

2. h{t) will serve as a Lagrange multiplier to define our problem as an initial condition 
problem. That is, we will fix h(t) such that the divergences from tadpole graphs 
contributing an inhomogeneous term to the equation of motion are cancelled and the 
initial value of the inflaton field is fixed for t < 0. Such a procedure will be implemented 
below when we obtain the equation of motion for the expectation value of the inflaton. 

Clearly, in De Sitter space-time the coupling to the curvature only serves to redefine the 
mass of the fields. For the inflaton, we just absorb this term in a redefinition of the mass, 
and set = without any loss of generality. 
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To study the time evolution of the expectation value of the inflaton field we invoke the 



tadpole method(see Ref. |19| and references therein) to obtain its equation of motion. This 
is implemented by first shifting $ by its expectation value in the non-equilibrium state as 
follows : 

$ ± (f, t) = <f>(t) + ^(x, t) ■ <t>(t) = ($ ± (f , t)) (3.3) 

and then enforcing the tadpole condition which is simply a consequence of ( p.3|) 

(iP ± (x,t))=0. (3.4) 

<p(t) being the spatially homogeneous inflaton zero mode that drives inflation. We also 
require that (a (x,t)) = to all orders, which means that the a field does not acquire an 
expectation value. After the shift of the inflaton field (|3.3j), the non-equilibrium action reads: 



S 



d 3 xdt \ C (iIj + ) + C (a + ) + ip + e 



3Ht 



3H<p — m%c 



(3.5) 



Jl 



(+ 



where £ represents the free theories for the ip and a fields. The tadpole condition (|3.4j ) is 
implemented by calculating the expectation value of the field operator ^ ± by inserting it into 
the non-equilibrium path integral and expanding in powers of g, about the free theory and 
setting the resulting expression to zero. In particular we evaluate the following expression 
to order g 2 in perturbation theory: 



(ip + (x,t)) = I P^P^le^W-^H) x ^ + (x,t) x 
oxp <; / d 3 x'dt' - 3Ht ' 



(3.6) 



) 2 + e 3Ht '^ + (a + ) 2 + e 3m '^- 



•)]} 







The equation of motion is obtained by independently setting the coefficients of (ijj + (x)tp + (x ; )) 
and '(x)ip~ (x 1 )) to zero. This is justified because these two correlators are independent 
functions. It is easily seen that the condition = yields exactly the same equations. 

This procedure gives rise to an equation of motion for the zero mode <fi(t) and is in fact 
equivalent to extremizing the one-loop effective action which incorporates the equilibrium 
boundary conditions at t < 0. 

It is worth noting that this scheme, referred to as amplitude expansion, also assumes that 
the amplitude of the zero-mode 4>(t) is "small" so that terms such as g<p(t)a 2 may be treated 
perturbatively. 

To one-loop order (0(g 2 )) we obtain the following form of the equation of motion. 
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m4> + ml<t>-Ua 2 ){t)-i^ f dTK(t,r) ( f>(T)- } ^- = 1 (3.7) 



2 ' w 4tt 

/oo 
^[G> 2 (r,t)-Gf(r,t)]e 3 ^, (3.8) 
o 

Spatial translational invariance guarantees that (a 2 ) is only time dependent. We can combine 
the independent contributions to the equation of motion into a single function: 

h(t) = ^ + 9 -(a 2 )(t). (3.9) 

If the a is conformally coupled to the curvature (a 2 ) represents a quadratically divergent piece 
which must be subtracted away by a renormalization of the "magnetic field" h. On the other 
hand if a is minimally coupled, one obtains the well known linearly growing contribution || 
after subtracting away the divergence leading to a modified inhomogeneous term in the zero 
mode equation of motion: 

7 , . gHH h(t) 

hit ) = £—- + -Li. 3.10 

W 8tt 2 2 V ; 

We will, however, show explicitly in the next section that this term does not enter in the 
calculation of the perturbative "decay rate" in this theory. 

We use h(t) as a Lagrange multiplier to enforce the constraint that for t < 0, <p(t) = 
<pi, <j)(t) = is a solution of the effective equation of motion for the expectation value. In 
this manner, we define an initial value problem with Cauchy data on a space-like surface for 
the dynamics in which the expectation value of the inflaton field is "released" from some 
initial value at time t = 0. For t > the equation of motion becomes 



0(t) + 3H(j)(t) + m%(j)(t) - I drK(t, r)0(r) = h(t) (3.1 1 : 



2 r t 







(j){t = 0) = (j>i , (j){t = 0) = (3.12) 



with the kernel given by 

We now proceed to study the solution to this equation of motion in the cases of conformal 
and minimal coupling. 



A. Conformally coupled massless a: 

For a conformally coupled massless scalar, the Green's functions are given by ( |2.27| , [2.28p 
with the Hankel functions as in (|2.26|) . 

The kernel has a particularly simple form in this case: 

^ r )=<*^ ffl . (3-13) 
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The /c-integral can be done by introducing a small convergence factor as follows 



sm(ak)dk = lim/m 
= lim 



e i(a+ie)k dk 



a 



*o a 2 + e 2 

Implementing the above prescription and absorbing a positive finite piece into e yields, 

H [ e H(t-r) _ X ] 



(3.14) 



K{t) T) l ~4 {e H ^))[(e H ^) - If + e 2 ] ' 
The equation of motion then becomes, 

g 2 E ft ^ [e H ^-l] 



(3.15) 



dr 



16tt 2 Jo (^-^[(e^*--) - l) 2 + e 



-Mr) = ~h{t). 



(3.16) 



A remarkable point to note, is that although time translation is not a symmetry in De 
Sitter space-time, nonetheless the term inside the integral is in the form of a convolution, and 
the equation, being linear in 0, can be solved via Laplace transforms. This is presumably 
due to the 0(4, 1) symmetry of De Sitter space. 

Notice that e serves as a short- distance regulator, and the e — > limit cannot be taken 
inside the integral because doing so results in a logarithmic short distance singularity in 
the equation of motion. This divergence is isolated by performing an integration by parts, 
keeping the e dependence, resulting in the equation of motion 



4>(t) + 3H<j)(t) + m|0(t) - 



g 2 H 

16tt 2 



H EE 



ln(l - e~ m ) 



-Ht 



I 



(r) 



+ ldr^at-r)\=h{t) 



(3.17) 
(3.18) 



£(t -r)= [ln(l - e' H(t - T ^) + e~ H ^ 
This expression makes it clear that the divergence can be absorbed in a mass renormalization, 



m% tR = m% + y|^(lne + l). 



(3.19) 



In what follows we will refer to m$ as the renormalized mass to avoid cluttering the notation. 
Taking the Laplace transform of the equations of motion, defining 



(j>{t)dt 



0( S ) = / , 

JO 

_ roc 
Z(S) = / ' 

Jo 

roc 

h(s) = / e~ st h(t)dt 
Jo 
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(3.20) 
(3.21) 
(3.22) 



and using the initial condition <p(t = 0) = 0, we find 



8 2 </>{a) - a<f>i + 3H{8(f>(s) - <f>i) 

dte~ st [\n(l - e- m ) + e~ m 



+ m% R 4>(s) + — —E(s)<j>(s) = h(s) 



Us) 



167T 2 

-si(s) 



(3.23) 
(3.24) 



B. The Gibbons-Hawking Temperature 

It is well known that fields in De Sitter space can be thought of, in some contexts, as 
being embedded in a thermal bath at a temperature H/2ir, the so-called Gibbons-Hawking 



temperature [24|. How do we see these effects in the context of our calculation? 



(2.23 


2.24 


. 2.27 


2.28 



,18| ) in ( |3.17| ) obey the periodicity condition in imaginary time 

27T2 N 

K(t-r) = K(t-r + —). 



(3.25) 



These kernels can be analytically continued to imaginary time t — r = — irj and expanded in 
a discrete series in terms of the Matsubara frequencies 



27m 1 H 

PH PH 27T 



(3.26) 



Th is recognized as the Gibbons- Hawking temperature in this space-time [23|. Then as a 



function of imaginary time ( 3.18Q can be written as 



PH n=—a 



We find 



S(uj n ) = -f3 H (— - S^j) , n > 1 



(3.27) 
(3.28) 

(3.29) 



^From this expression, E(s) can be obtained at once: 

1 1 , 



t(a) = -ai(8)= y £H( 



nH s + nH s + H 



(3.30) 



This Laplace transform has simple poles at minus the Matsubara frequencies corresponding 
to the Hawking temperature. 

Finally we obtain the solution for the Laplace transform 
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- = »(. + 3H) + fc(.) (331) 
S 2 + 3i/.s + m| + 

This expression clearly shows that the information on the "decay rate" is obtained from 
the imaginary part of the poles of the denominator and that the "effective magnetic field" 
h(t) whose Laplace transform h(s) is explicit in eq. (|3.31 ) is not relevant to understand this 



"decay rate". In fact as we will show shortly, the "decay rate" is determined only by the 
"self-energy" S(s). Since our goal is to understand the "decay rate" we set h(s) = without 
loss of generality. 

The denominator in the above expression is the inverse propagator for the massive scalar, 
evaluated at zero spatial momentum and S(s) is the one-loop correction to the self-energy. 
However, from ( |3.3Up we see that S(s) is analytic in the s-plane with simple poles at s — 
—nH, with n an integer^ 1). This is a surprising result, because we expect that, as in 
Minkowski space, there would be a cut in the s-plane, which since the particles in the loop 
are massless, would run along the entire imaginary axis. Obviously, this is not what happens 
- the analytic structure is very different in De Sitter space-time from that in Minkowski 
space-time. In fact we find that <f>(s) has two simple poles, perturbed from their original 
positions by 0(g 2 ) corrections. We will understand this result by taking the flat-space limit 
in a later section. 

The real time dynamics of the zero mode of the inflaton is now found by inverting the 
Laplace transform by complex integration: 

m = / e st <j>(s) — (3.32) 

where the Bromwich contour is taken to the right of all the singularities. 
The poles of <fi(s) can be obtained by setting, 



s 2 + 3#s + m| + -|— S(s) = 0. (3.33) 

We analyze the 4m| > 9H 2 and 4m| < 9H 2 cases separately. In the first case, with 
4m| > 9H 2 , we have that in the absence of interactions the inflaton undergoes damped 



2 



oscillations red-shifted by the expansion. At zeroth order the poles (|3.33|) are at 



± -3H ,yj4m%-QIP 



(3.34) 



We absorb the real part of the self energy on shell in a further (finite) redefinition of the 
renormalized mass. Since S(s) is a meromorphic function of s, i?e[E(so )] = -Re[S(sg )]• Then 
to this order the "pole mass" is given by 



m f = m l + g2Re }^ 0n . (3.35) 
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Obviously this mass is independent of the renormalization scheme. We also define the sub- 
tracted self energy S*(s) = S(s) — i?e[S(so")], and at this stage it is convenient to introduce 
the "effective mass" 



Ml 



*2 

m$ - 



9H 2 



(3.36) 



in terms of which, to order g 2 , we find 

_u -3H 



2 

-3if 



(3.37) 



The Bromwich countour for the inverse Laplace transform can now be deformed by 
wrapping around the imaginary axis and picking up the poles . The resulting expression 
has the Breit-Wigner form of a sharp resonance centered at the "pole mass" that is very 
narrow in the weak coupling limit, leading to the time evolution of the inflaton given by 



<p{t) = Z4>ie- 3Ht/2 e- rt/2 cos(M $ t + a) 



(3.38) 



where, 



1 + 



9H 2 



2g 2 B' 



12g 2 HB 18g 2 H 2 B' 



4Mg 32vr 2 M $ 128vr 2 M| 



128tt 2 M| 



a = arctan 



9' 



3H 



2M$ ' 16vr 2L 2M $ ' 8M| 



A' 9H 2 A' 
+ 



18H 2 B 
16M| ^ 



(3.39) 



B = Im 



A' = Re 



B' = Im 



S*(s 

ds 

dt*(s] 
ds 



s==fl-+iM <s , 



-3H 
2 



The most interesting of these quantities is the damping rate, given by 

g 2 Imt*(s+) 



r = 



16tt 2 M$ 



(3.40) 



The imaginary part of the self-energy in the above expression, can be evaluated using stan- 
dard sum formulas (these sums also appear finite temperature field theory H25,26|) and is 
found to be 
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This expression for T is a monotonically increasing function of H. When H = 0, it matches 
with the flat-space decay rate 



9 2 

^Minkowski TTTi " (3.42) 

We thus obtain the result that the inflaton decay proceeds more rapidly in a De Sitter 
background than in Minkowski space-time. 

The rate ( |3.41| ) can be written in a more illuminating manner as 

r _ ^coth(te) g* [1 + 2n b {«, )] 

327rM 4 32ttM$ 1 ' ' 

u = -is+ ; n b (uj ) = e0HUQ _ - (3.44) 

This is a remarkable result - the decay rate is almost the same as that in Minkowski space, but 
in a thermal bath at the Hawking temperature |19| , [27|] . Habib |15| has also found an intriguing 
relationship with the Hawking temperature in the probability distribution functional in his 
studies of stochastic inflation. 

If we now look at the 4m| < 9H 2 case, we see that the inflaton ceases to propagate, i.e. 
its Compton wavelength approaches the horizon size and there is no oscillatory behavior in 
the classical evolution of the zero mode: 

= (^ - l)- 1 /20. e -3^/2 ginh (| M$ | t + ^ ( 3 _ 45) 

4m| 

tanh 3 <J 



3H 



The one loop contribution is the same as in the previous case but now the poles of <f)(s) lie 
on the real axis: 



, -3H J9H 2 -4ml 

5? = ±- . 

2 2 



(3.46) 



s 



± _ -3H J9H* - 4m% - 4£,£(s±) 



We define 

S(4) = C±D (3.47) 
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and absorb C into a finite mass renormalization, 



a 2 C 

mf = m% + y— ^. (3.48) 
Inverting the transform we obtain, 

<j>(t) = Z^e' mt l 2 e- Vt l 2 sinh A + /3), (3.49) 
with Z being the (finite) wave function renormalization i.e. the residue at the poles, and 



K = ^9H 2 - Amf (3.50) 



g 2 tan(W9 - ^f-) ™* On?* 

r = y 2 V H2 for ^ >jfJ> f!^. (3.51) 



1Qtt^/9H 2 - Amf V2 



We also see that the numerator of the damping rate in eq.(3.46) diverges at H = ^| > ^y 1 -, 
which indicates the breakdown of perturbation theory. 

The decay rate may again be written in the a form that makes explicit the effect of 
the Hawking temperature and the "stimulated decay" with the Bose-Einstein distribution 
function at the Hawking temperature. 

' ^ 1 + 2 ^°)] ;wb = -«f. (3.52) 



m% 2 



We can also study how the "decay" of the inflaton into lighter scalars modifies the "slow- 
roll" evolution such as would occur in new inflationary models. Since we are assuming a 
small field amplitude for our pertubation theory to make sense, we cannot treat chaotic 
inflationary scenarios in the same way. 

For this we now set m| < and study the situation in which the inflaton is "rolling" 
from the top of the potential hill, during the stage of quasi-exponential expansion. In this 
situation the amplitude expansion is valid when the zero mode is close to the origin and at 
early times. 

In this case it is convenient to write m% = — /z 2 . Now the poles are on the real axis at 
the positions: 



- 4 =p ' - ,3,4) 



V 2 / 



, -3H , U3H 

Sn = — — ± 
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For g = the pole that is on the positive real axis, Sq , is the one that dominates the 
evolution of the inflaton down the potential hill (the growing mode). Since Sq > we find 
that 

oo + 

^» +) = ?^(4T^) >0 - (3 ' 56) 

and we conclude that, to this order, the pole on the positive real axis is shifted towards 
the origin. Therefore the rate of growth of the growing mode is diminished by the decay 
into lighter scalars and the "rolling" is slowed down. This is physically reasonable, since 
the "decay" results in a transfer of energy from the inflaton to the massless scalars and the 
rolling of the field down the potential hill is therefore slowed down by this decay process. 



C. Minimally coupled massless a: 



Now we turn to the £ CT = case. The Green's functions are given by ( |2.27|J2.28|) with the 



Hankel functions fl2.25p . There are many features in common with the previous case: the 
kernel is translationally invariant in time and periodic in imaginary time with periodicity 
(3h- Thus it can again be expanded in terms of Matsubara frequencies and the Laplace 
transform carried out in a straightforward manner. However, for minimal coupling we find 
a new logarithmic infrared divergence in the kernel along with the ultraviolet logarithmic 
singularity of the conformally coupled case. Therefore, the /c-integral in the kernel must be 
performed with an infrared cutoff /i. The logarithmic ultraviolet divergence is handled in as 
before, and after this subtraction we find the self-energy to be: 

£ CT =l/6 



°°1 1 q AJ-f 00 1 1 1 

= E - — —) - -"vE -(— -JT ' — ~s) (3-57) 

j nH s + nH s + H 6 ^ n s + nH s + SH + nH 

AH 1 1 N AH 2 r l 1 



3 K s + H s + 2H' 3 L s 2 (s + 3H) 21 ' 

We notice that the self-energy for the conformally coupled case is contained in the above 
expression. It is also interesting to note that the imaginary part of the self-energy, on-shell, 
receives no contribution from the infrared divergent piece, which only contributes to a further 
renormalization of the "pole mass" m*£. We now focus on the case 4m| — 9H 2 > 0. 

The expression for the decay rate is obtained by evaluating the imaginary part of the 
self-energy at the pole Sq (see eq. ( |3.40| ) and found to be 




r=^[^7^|l + ^]+^]- (3-58) 



18 



with ( |3.44D . The qualitative behaviour of the decay rate as a function of H is similar to that 
in the conformally coupled case. We see that, r(£ CT = 0) > r(£ CT = 1/6) > r (Minkowski). 



D. Analytic structure of the self-energy 

As mentioned earlier, the analytic structure of the self-energy in the De Sitter background 
is drastically different from what one would expect in flat space-time. We find that the self- 
energy for the case where $ is unstable does not display cuts in the s-plane. Though the 
reason for this strange behaviour is not clear at present, we can try to understand it by 
taking the flat-space limit i.e. the H — > limit. In this limit, it suffices to look at E(s) for 
£ CT = 1/6. It is easily shown that the additional terms in (3.25) are vanishingly small in the 
H — > limit. We thus obtain, 



lim S(s) = limfV HI 



1 



h^o h^o 1 ^ nH s + nH s + H 

roo I 1 



H X S + X 



) dx — 1 



(3.59) 
(3.60) 



= ln(s/m) + ln(m/H) — 1 
where m is a mass scale that serves as an infrared regulator. Hence, 



(3.61) 



s 2 + (m% + £s Hm/H) - 1) + ^ \n(s/m) 

The above expression has a cut which can be chosen to run from to — oo, so that it is to 
the left of the Bromwich contour. The discrete singularities of the self-energy at s — —nH 
have merged into a continuum in the small H limit to give a cut along the negative real 
axis in the s-plane. In addition, there are two simple poles in the first Riemann sheet at 
±im$ + 0(g 2 ). Inverting the transform yields, 

<f)(t) = Z4>ie~ rt/2 cos(m<s,t + a) 

— (pi / dui 5 5 — . (3.62) 

We see here, a Breit-Wigner form plus contributions from across the cut. 

A flat-space analysis with m a ^ 0, carried out in a previous work reveals a very 



different picture. In particular, we find two cuts extending from s = ±2im a to ±ioo and two 
simple poles which move off into the second Riemann sheet above the two-particle threshold 
when m$ > 2m (T . However, upon inverting the transform and taking m a — > 0, we obtain 

^16Wo W (^_ m |__^i n(a; / m )) 2 + (£)2- 
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A deformation of the contour of integration shows that eqs. (|3.62] ) and ( |3.63| ) are in fact 
identical and hence we have consistency (see the Appendix). Although the limits H — ► 
and rricr — > do not commute insofar as the analytic structure in the s-plane is concerned, 
the time evolution obtained from the inverse Laplace transform is unambiguously the same. 



IV. THE 0(2) MODEL AND GOLDSTONE BOSONS 

While the inflaton is typically taken to be a singlet of any gauge or global symmetry, there 
is no reason that it could not transform non-trivially under a continuous global symmetry. 
If this occurs, and the symmetry is broken spontaneously, the interesting possibility arises 
of dissipation of energy into Goldstone bosons. Even if the inflaton is a singlet, one could 
ask about the evolution of fields that belong to multiplets of a global continous symmetry, 
during inflation. Our primary motivation, however, is to study dissipative processes of the 
inflaton via the decay into Goldstone bosons in De Sitter space-time as new non-equilibrium 
mechanisms. 

The effect may be studied in the 0(2) linear sigma model, in which the inflaton is part 
of the 0(2) doublet, with Lagrangian density 

C = i[e 3ffi $ 2 - e m (V$) 2 + e 3Ht n 2 - e m (Vvr) 2 ] (4.1) 
+I e 3H*( M 2 _ i2£# 2 )($ 2 + n 2 ) - e 3m ^($ 2 + 7T 2 ) 2 + h(t)$. 

Assuming that /i 2 > 12^ H 2 , in the symmetry broken phase the field $ acquires a VEV 



(*> = f" 2 -^' . (4.2) 
and the Goldstone bosons are left with no mass terms, and hence are massless, minimally 



coupled fields in De Sitter space-time p8| ; this is a fairly well known result. 
We now write 



t) = y 6 ^ 2 l ^ R2 ) + m + t) (4.3) 

and use the tadpole method to impose: 

(#M)> = 0, (4.4) 
(n(x,t)) = 0. 

Using eq.(4.3) it is easy to see that couplings of the form tpn 2 , are automatically induced in 

the Lagrangian. Carrying out the same analysis as before and keeping only terms up to first 
order in 4>(t), we find the following equation of motion 
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+ 3#0 + m 2 (t)0 - f dr[K^t, r) + ]-K„{t, r)]0(r) = h(t), (4.5) 

/oo 
k 2 dk\Gf (r,t) - Gf(T,t)}e 3HT . 
o 

Here m 2 (t) is given by the tree- level mass term 2/z 2 plus contributions proportional to 
the tadpole (7r 2 )(t). Since the Goldstone boson is always minimally coupled to the curvature 
(7r 2 )(t) will grow linearly with time and result in a non-trivial time- dependence for the mass 
of the massive mode <fi. In this case the methods of the earlier section are not applicable, 
and the definition of a "decay rate" will be plagued by time-dependent ambiguities. 

Such a time-dependence will, however, be absent in a De Sitter space that exists forever, in 
which case the contribution from the tadpole will simply be a divergence that is renormalised 
through a redefinition of the mass. In such a situation, the contribution to the decay rate 
from the massive mode would be subdominant due to the kinematics and the Goldstone 
modes would provide the dominant contribution to the decay rate. Thus, neglecting K^t, r) 
the equation of motion has exactly the same form as eq.(3.5). The subsequent analysis and 
results would be identical to those obtained above in the section on minimally coupled fields. 



V. THE LANGEVIN EQUATION 

Our results allow us to make contact with the issue of decoherence and the stochastic 
description of inflationary cosmology [|T^ -|15ll. 

Decoherence is a fundamental aspect of dissipational dynamics and in the description of 
non-equilibrium processes in the early universe [p^ -[T5|]. The relationship between fluctuation 
and dissipation as well as a stochastic approach to the dynamics of inflation can be explored 
by means of the Langevin equation. 

This section is devoted to obtaining the corresponding Langevin equation for the non- 
equilibrium expectation value of the inflaton field (zero mode) in the one loop approximation 
within the model addressed in this article. 

The first step in deriving a Langevin equation is to determine the 'system' and 'bath' 
variables, and subsequently integrate out the 'bath' variables to obtain an influence functional 
for the 'system' degrees of freedom [|T3|-[15]. 



We first separate out the expectation value and impose the tadpole conditions as follows: 

d>± = 0± + ^ ; ($ ± (f,t))=0 ± (t) ; (a ± (x,t)) = 0. (5.1) 

The non-equilibrium effective action is defined in terms of the Lagrangian given by equa- 
tion(|3.2|) as, 

i{Sol<f>+]-S [<f>-]}J{S iip+]-S li>-]} 



J Vi) + Vi)-Va + Va- e i{So[,t>+] - s ^ ]} e 



e i{S [<T+]-S [a-]} e i{S in t[i>+,cl>+,a+]-S int [f-,cl>-,a-]} 
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where So represents the action for free fields, and Si nt contains all the remaining parts of the 
action. The effective action for the zero modes is obtained by tracing out all the degrees of 
freedom corresponding to the non-zero modes (or the 'bath' variables), in a consistent loop 
expansion. It is useful to introduce the center of mass (</>(t)) and relative (R(t)) coordinates 



as 



= <f>(t) ± 



R(t) 



(5.2) 



Remembering the definitions of the Green's functions in eqs.( |2.29| - |2.31| ), we expand 
exp(iS int ) up to order g 2 and consistently impose the tadpole condition to obtain the ef- 
fective action per unit volume (because the expectation value is taken to be translationally 
invariant), 



S, 



eff 



n 



dt [Co 

2 



167T 3 

ig 2 



CO 



dt 

— CO J —CO 

oo rco 



dt 



32tt 3 

+ h(t)R{t) + terms independent of 



dt'Jd 3 ka 3 (t)a 3 (t')R(t)<f)(t') [G> 2 (t,t')-Gf(t,t') 
dt' [ d 3 ka 3 (t)a 3 (t')R(t)R(t') [G> 2 (t, t') + G?(t, f 



(5.3) 



As before h(t) results from subtracting the (^-independent contribution of the tadpole 
from h(t). 

Now, it is not difficult to verify that the non-equilibrium Green's functions given by 
eqn.( [2.16| , p,19[ ) and obeying the boundary conditions Q2.22j ) have the property, 

G>(t,t') = [G<(t,t')r. (5.4) 

This property guarantees that the second term in the effective action is real, while the third 
term is pure imaginary. The imaginary, non-local, acausal part of the effective action gives 
a contribution to the path-integral that may be written in terms of a stochastic field as 

exp[-^jdtj dt'R(t))C(t, t')R{t')](x J V£P[i] exp[z j di£(t)R(t)\ (5.5) 

m = e M-\J dt J dtf^t)Kr\t, oecoi 



with 

JCM = -^Jd 3 ka 3 (t)a 3 (t'){G> 2 (t,t>) + G< 2 (t,t')] (5.6) 
The non-equilibrium path integral now becomes 
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Zoc J V(f>VKD£V[£\ expi[S re// (0, R) + J dt£(t)R(t)] (5.7) 

with S re ff[4>, R) being the real part of the effective action and with V[£] the gaussian prob- 
ability distribution for the stochastic noise variable. 



The Langevin equation is obtained via the saddle point condition [13] 



6Sreff \ R =o = m (5i 



5R(t) 



leading to 



+ 3i j0 + m |0 + ,^!_ f dt ' [d 3 ka 3 (t'){G> 2 (t,t') - Gf (M')]<KO -Ht) = (5.9) 

107T J — oo J a ( It J 

where the stochastic noise variable has gaussian, but colored correlations: 

« »= 0; « £(*)£(*') »= K(t, t') (5.10) 

The double brackets stand for averages with respect to the gaussian probability distribution 
V[£\. The noise is coloured and additive. Due to the properties of the Green's functions (5.4), 
the noise kernel (which contributes to the imaginary part of the effective action) is simply 
the real part of G^ 2 (t,t'). The dissipative kernel (which gives rise to the non-local term 
in the Langevin equation) on the other hand is given by the imaginary part of G^ 2 (t,t ; ). 
Consequently, the fluctuation — dissipation theorem is revealed in the Hilbert transform 
relationship between the real and imaginary parts of the analytic function G^ 2 {t,t'). The 
equation of motion ( |3.11| ) is now recognized to result from ( |5.9| ) by taking the average over 



the noise. 

One could now obtain the corresponding Fokker-Planck equation that describes the evo- 



lution of the probability distribution function for <f> [|12l , p9 • 

We must emphasize here, that all the results discussed so far in this section, are inde- 
pendent of the specifics of the FRW background spacetime, the masses of the fields, and the 
temperature of the initial thermal state. The importance of the Langevin equation resides 
at the fundamental level in that it provides a direct link between fluctuation and dissipation 
including all the memory effects and multiplicative aspects of the noise correlation functions. 

In particular in stochastic inflationary models it is typically assumed |12| that the noise 



term is gaussian and white (uncorrelated) . This simplified stochastic description, in terms 
of gaussian white noise, leads to a scale invariant spectrum of scalar density perturbations 
fl~2|l . Although this description is rather compelling, within the approximations made in 
our analysis we see that for a a field with an arbitrary mass and arbitrary couplings to 
the curvature there is no regime in which the correlations of the noise term ( |5.10| ) can be 
described by a Markovian delta function in time |T3|]. One could speculate that some other 
couplings or higher order effects, or perhaps some peculiar initial states could lead to a 
gaussian white noise, but then one concludes that gaussian white noise correlations are by 
no means a generic feature of the microscopic field theory. 
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We can speculate that our result of gaussian but correlated noise could have implica- 
tions for stochastic inflation. In particular it may lead to departures from a scale invariant 
(Harrison-Zel'dovich) spectrum of primordial scalar density perturbations, that can now be 
calculated within a particular microscopic model using the non-equilibrium field-theoretical 
tools described in this work. Of course this requires further and deeper study, which is 
beyond the scope of this article. 

VI. CONCLUSIONS AND FURTHER QUESTIONS 

In this article we have studied the decay of inflatons into massless scalars in De Sitter 
space. The motivation was to understand the non-equilibrium mechanisms of inflaton relax- 
ation during the stage of quasi-inflationary expansion in early universe cosmology, in old, 
new and stochastic inflationary scenarios. Chaotic inflation models cannot be treated within 
our approximation scheme since the field will typically start at values that are large enough 
so that the amplitude expansion is not valid. Having recognized the inherent difficulties in 
defining a decay rate in a time dependent background, we have given a practical definition 
that requires understanding the real-time evolution of the inflaton interacting with lighter 
fields. This led us to a model in which an inflaton interacts with massless scalars via a 
trilinear coupling, allowing for both minimal and conformal coupling to the curvature. 

We obtained a decay rate at one loop order that displays a remarkable property. With 
some minor modifications it can be interpreted as the stimulated decay of the inflaton in a 
thermal bath at the Hawking temperature. This decay rate is larger than that of Minkowski 
space because of the Bose-enhancement factors associated with the Hawking temperature. 
The decay rate of the minimally coupled case is larger than that of the conformal case which 
in turn is larger than the Minkowski rate. We have shown explicitly that in the case of new 
inflation, the dissipation of the inflaton energy associated with the decay into conformally 
coupled massless scalars slows down further the rolling of the inflaton down the potential 
hill. We have also studied decay into Goldstone bosons. 

To establish contact with the stochastic inflationary scenario, we derived the Langevin 
equation for the coarse-grained expectation value of the inflaton field to one loop order. 
We find that this stochastic equation has a gaussian but correlated (colored) noise. The 
two point correlation function of the noise and the dissipative kernel fulfill a generalized 
fluctuation-dissipation relation. 

There are several potentially relevant implications of our results for old, new and stochas- 
tic inflation. In the case of old inflation we see that dissipative processes in the metastable 
maximum can contribute substantially to the "equilibration" of the inflaton oscillations, and 
more so because of the enhanced stimulated decay for large Hubble constant. In the case of 
new inflation, we have seen that these dissipative effects help slow the rolling of the inflaton 
field down the potential hill, possibly extending the stage of exponential expansion. 

Within the context of stochastic inflation, our results point out to the possibility of 
incorporating deviations from a scale invariant spectrum of primordial scalar density pertur- 
bations by the noise correlations, which manifest the underlying microscopic correlations of 
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the field theory. This is a possibility that is worth exploring further. 
Our results also point to further interesting questions: 

1. Is it possible to understand at a more fundamental level the connection between the 
decay rate and the Hawking temperature of a bath in equilibrium! In particular, is 
this connection maintained at higher orders? 

2. Is it possible to relate the "decay rate" to the rate of particle production via the 
interaction? Such a relation in Minkowski space-time is a consequence of the existence 
of a spectral representation for the self-energy but such a representation is not available 
in a time dependent background. 

Answers to these questions will undoubtedly offer a much needed deeper understanding 
of non-equilibrium processes in inflationary cosmology. 



ACKNOWLEDGMENTS 

D.B would like to thank H. J. de Vega for useful references and illuminating discussions, 
and NSF for support through grant award No: PHY-9302534. D. B. and S. P. K. would like 
to thank Salman Habib for illuminating conversations and comments. 

S. Prem Kumar and R.Holman were partially supported by the U.S. Dept. of Energy 
under Contract DE-FG02-91-ER40682. 



VII. APPENDIX 



In this appendix, we show explicitly, that ( |3.62j ) and Q3.63Q are in fact identical expres- 
sions. The behaviour of 4>(t) in Minkowski (given by (3.57)) was obtained in reference |19 | 



where m a was allowed to have non-zero values. Now, the results from this reference must 
obviously coincide with the flat-space limits (H — > 0) derived in this paper. To see this, 
consider eqn. (|3.63|) : 



g 2 f°° u cos out 



-ReiT -iduu(uje iut ){- 

71 J L 



l-^H^M+i^ 



uj* — m 



The branch cut due to the logarithm in the denominator runs from to +oo, and the contour 
of integration may be chosen to run from to +oo on the first sheet. It is easy to see that 
the integrand has 4 simple poles, one in each quadrant. In particular, the pole in the first 
quadrant is at 
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g 2 g 2 

uj 1 = m<£ + —— ln(ra$/m)+i— . (7.1) 

327r^m$ 647rm$ 

Furthermore, the exponential is well behaved at infinity in the first quadrant. We can 
therefore deform the contour of integration so that it runs /,from to +200, picking up the 
residue from one simple pole. Writing oo = iz, and noticing that the logarithm picks up an 
imaginary part m /2 we now have, 



(f)(t) = -RelT 'dz(ize- zt )[ 

IX J 



1 



-z 2 — m% — -^\n(z/m) - i-^ 

iue 



167T 2 V / ) 167T 

+ Res[- 



-z 2 -ml- ^\n(z/m) u 2 - m% - ^\n(u/m) - i£ 



-y2 foo Z6~ Z ^ 

dz 



: r 9 f a 



+ 2ixRe\ lim \ ^ - 

_ m | _ ^In^/m) - i£ JU 

Absorbing the real part of the self-energy on-shell as an additional mass renormalization, 
and going through the subsequent algebra we obtain our result, 

q 2 ^ 
<f>(t) = 0(0) (1 + 327r 2 m 2 ) e M ™* cos ( m ^) 

~ ( t >i ~Tn — 9 / ^ 



1Qtt 2 J o (z 2 + m| + ^_\n(z/m)) 2 + 



16vr' 
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